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Abstract 

In the interaction of two electric quadrupoles, there is at short 
range a contact term proportional to the second derivative of a delta 
function. This contact term contributes to the hyperfine splitting of 
bound states of two particles with spin one or higher - for example 
the bound state of 0,^ and a nucleus with spin one. The contact 
hyperfine splitting occurs in states with orbital angular momentum 
unity (p-wave), in contrast to the Fermi contact interaction which is 
in S'-states. We find that these contact splittings will be observable 
with 0,~ atoms and help measure the quadrupole moment and charge 
radius of the hyper on. 
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1 Introduction 



It is well known in the context of atomic liyperfine interactions that the 
magnetic dipole-dipole interaction has a contact term, called the Fermi con- 
tact interaction Ij. A well known example is the splitting of F = from 
the F = 1 states of the H atom, responsible for the microwave radiation 
of interstellar Hydrogen. It is natural to ask whether such contact interac- 
tions occur in higher multipoles, for example in the interaction of two electric 
quadrupoles. Only particles with angular momentum J = 1 or higher can 
have quadrupoles. We consider here as an example the bound state of an Q~ 
hyperon ( J=3/2 ) to a nucleus of angular momentum one (say NitrogenlA), 
though there are other possibilities like a bound state of an anti-deuteron to 
a nucleus, etc. We find in addition to the well known long-range electrostatic 
interaction of two quadrupoles (~ there are short range contributions 

some of which are proportional to the second derivative of a delta func- 
tion such as A6{r). These contact terms give rise to splitting in the P- wave 
(L = 1) states of the bound system of an Q~ with a nucleus, or other systems 
with two quadrupoles. As far as we know this higher contact interaction has 
never been described, although it is a consequence of ordinary electromag- 
netism and quantum mechanics. We find that this splitting may be useful to 
determine the electric quadrupole moment of the Q~ , which is not yet known 
experimentally. 

In an ordinary atom, an electron is bound to the nucleus. The electron 
has spin, and associated with the spin there is a magnetic moment. This 
magnetic moment interacts with the magnetic fields in the atom and leads to 
splittings in the energy levels which are called fine structure- if the nucleus 
itself has no magnetic moment, for example a spinless nucleus such as Helium 
four. However if the nucleus itself has spin and magnetic moment, the nu- 
clear magnetic moment (which is much smaller than the electronic magnetic 
moment) gives rise to further splittings which are called hyperfine structure. 
Other contributions to atomic hyperfine structure arise from nuclear electric 
quadrupole moment ( which is sensitive to electric field gradients at the nu- 
cleus). Higher nuclear moments tend to produce smaller effects. In an fl~ 
atom where the Q~ is bound to a nucleus, there are interesting effects com- 
ing from the magnetic moment of the (which is known experimentally) 
and the electric quadrupole moment (which is not known experimentally but 
estimated theoretically). The fine structure in an Q~ bound to a spinless 
nucleus was discussed by Goldhaber and Sternheimer j2]. In the scenario 
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of Ref. [2] the electric quadrupole dominates over magnetic contribution to 
the fine structure because the electric quadrupole is assumed to be large. 
However later estimates of the Q~ quadrupole, using QCD inspired quark 
models [Hj are about a hundred times smaller than assumed in reference 
This leads to a completely different scenario, in which the fine structure is 
dominated by magnetic interactions, and the effect of the electric quadrupole 
is to be found in hyperfine interactions. This is the scenario we discuss here 
by taking a nucleus which has an electric quadrupole and a (small) magnetic 
dipole to bind the Q~. The hyperfine effects will now be smaller than fine 
structure but still large enough to be observed. As far as we know these 
hyperfine effects have not been discussed in the literature. The most inter- 
esting of these effects is in the interaction of two electric quadrupoles, where 
there is a contact term in the / = 1 (p-wave) states. This contact term is 
discussed in Section 2 and Appendix 1,2. Then in Sections 3, 4 we discuss 
an application to fl~ atoms. We conclude in Section 5 with a summary of 
the results. 



2 Contact interaction of quadrupoles. 

The long range interaction of two quadrupoles Q^^^ and Q*-^-* may be written 
in terms of Cartesian components in the form 

rltensor - ^g^^ij ^ j{ ^5 ^kl ) " 

= i^[2gSj)gSf - 20QS;)gS?n,n, + 35QS])Qi?n.n,n,n,] (2.1) 

In equation (2.1) repeated indices are summed; in the first term the two 
quadrupoles are fully contracted and in the next terms rii = Xi/r, etc. The 
interaction of two quadrupoles decreases very rapidly with distance- like r~^, 
as follows from dimensions- much faster than the fine structure coupling, 
which decreases only like r~^. As a result the quadrupolar coupling is more 
important for smaller orbits, at low quantum numbers. Therefore we will 
focus on low lying P-states of the atom. It can be checked that the tensor 
quadrupolar interaction (eq. 2.1) is ambiguous in P-states of the atom; the 
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angular averaging vanishes because the angular dependence of the interac- 
tion (when expressed in spherical harmonics) involves Y4^rn{(^ , 4>) , where the 
angles 6', are those of the vector connecting the two quadrupoles ,4,. The 
expectation value of a spherical harmonic Y^^m vanishes in a P-wave state. 
The radial integral on the other hand diverges. This ambiguity is similar to 
the case of the dipole-dipole interaction in S-waves, where the contact inter- 
action solves the ambiguity. Similarly here, instead of the tensor interaction 
of two quadrupoles we have a contact term, involving the derivative of a delta 
function, which generalizes the contact interaction where there is only a delta 
function. The quadrupolar contact coupling is a scalar interaction coupling 
the two quadrupole tensors and second derivatives of a delta function (as 
required by dimensional considerations) to a scalar. The precise formula for 
the quadrupolar contact term is given below 

Heont = -^[^(V.V, - lv%,)6^'\r) + ^6,,V'6^'\r)]Q^^Q\l\ (2.2) 

Equation (2.2) is derived in the Appendix 1. In this formula Q*^^-*'*^^^ are 
the quadrupole moments of the hyperon and the nucleus, and = VjVj 
is the Laplace operator. This contact term gives a finite contribution in a 
P-wave state where the two spatial derivatives are taken up by the wave 
function squared. 



3 Structure of Q atom 



Consider an interacting with nucleus with charge Ze. In the first approx- 
imation this is Coulomb system with spectrum 

where n, / are principal and angular momentum quantum numbers, m is the 
reduced mass 

m 



[1 + mn/M) 

ruQ is Q mass and M is the mass of the nucleus. The corresponding Rydberg 
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constant for heavy nucleus is 



= rJ!^ = AAkeV. (3.2) 

The interaction of charge with magnetic and quadrupole moments and of 
the multipole moments with each other is smaller than Coulomb interaction 
of charges and can be considered as perturbation to Coulomb levels eq.(3.1). 
In the case of ordinary atoms perturbations split Coulomb levels into two 
main structures - fine structure (due to interaction of magnetic moment of 
electron with angular momentum of atom) and "then" fine structure levels 
into hyper-fine structure ( due to interaction of electron magnetic moment 
with nuclear magnetic moment). In the case of Q~ atom the splitting of 
the Coulomb levels can be more complicated because Q~ has spin S = 3/2 
and magnetic and quadrupole interactions can have rather different scale of 
strength. It is not a big deal to write equations for splitting for any possible 
case. In this paper we restrict ourselves to the case most similar to the 
ordinary atoms, i.e. we suppose that there is fine structure due to interaction 
of magnetic moment of Q~ with orbital angular momentum L. That means 
in particular that magnetic moment of nucleus is much smaller than the 
magnetic moment of and that corresponding interaction of quadrupole 
moments with L is also much smaller. 

For fine structure levels we have one "good" quantum number-total an- 
gular momentum J that is the sum of orbital angular momentum L and spin 

s of n--. 

J = L + S. 

Each Coulomb level E^^i splits in general into four levels with different 
total angular momentum j = 1 + 3/2, j = 1 + 1/2, j = l — l/2,j = 1 — 3/2. The 
magnetic interaction for Q~ atoms is similar to the case of ordinary atoms 

^-.--^(^-1)L-S, (3.3) 

where factor g relates magnetic moment with spin /xn = g^- 

The average value of has to be calculated with unperturbed Coulomb 
wave functions 

< — >= . (3.4) 

r3 n3/(/ + l/2)(/ + l) ^ ' 

Thus the fine splitting is equal to 
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AE{n,l,j)- ^3 m[g-l) 2/(2/ + 1)(/ + 1) ' ^^'^^ 

Magnetic moment of Q~ is known experimentally fiQ = —2.02 ± 0.05 /Uat 
Thus substituting these factors into Vmag we get for fine splitting of p-wave 
levels into three states with total angular moment j = 5/2, 3/2, 1/2 : 

^4 23 
AE{n, = Y^[J(J + 1) - ^^-^^ (^-^^ 

This formula gives the scale and exact form of the fine splitting. 

We are interested in the hyperfine splitting of these levels by interaction of 
Q quadrupole moment with nucleus. In this case only total angular moment 
of the system F = J + I is a good quantum number. (Here I is a nuclear 
spin). The splitting due to long range interaction was considered in ref. 
But interaction of quadrupole moment with charged system has also local 
terms proportional to 5-function and to it derivatives. Thus there should be 
additional hyperfine structure for states with small angular momentum, i.e. 
additional small shifts of s- and p- states. We are interested in hyperfine 
structure due to these contact terms. 

There are well-known local terms, such as Fermi contact interaction of 
two magnetic moments /Xi 2 

f/(x) = -^(/irVf^)5^'nr). (3.7) 

Fermi term gives contribution into s-wave only and it does not compete with 
hyper-fine structure due to quadrupoles in p-waves. 

Interaction of charge ei,2 with e.m. radius r2,i and radius-radius interac- 
tion are also local 

U{x) = -e,e2[krl + rl)6^'\v) + LrlrlV'6^'\r)]. (3.8) 
o 3o 

First term in eq.(3.8) contributes into s-wave states shift, the second one 
into s- and p-wave states. But both of this terms do not depend on the spin. 
Thus they contribute into shift but not into splitting of states with different 
J and/or F. 

Contact terms from eqs. (3.7) and (3.8) are well known and their deriva- 
tion can be found in text-books on Classical Electrodynamics ( see, e.g., ref. 
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0). Similar local terms exist also in the interaction of quadrupole moments, 

which we cannot find in the hterature. 

There are two such terms ( for derivation see Appendix 1): 
Interaction of quadrupole moment Qij with e.m. radius r 

Vix) = -^[eirfQlf + e2rlQ[f]VNjS^'\rl (3.9) 
and quadrupole-quadrupole interaction 

V{x) = -^Qg^QS?[V,V, - lv^5,,]5(^)(r). (3.10) 

These contact terms contribute to sphtting of s- and p-waves states. 

Notice that long-range interaction of two quadrupoles does not con- 
tribute into splitting of s- and p-states. Any matrix element of this ci-wave 
interaction vanishes identically for s- and p-states. 

We are interested in determination of the value of Q quadrupole moment 
from the spectrum of Q atom. Thus we concentrate on hyperfine structure 
of p-states. 

There are three sources of splitting for p-states as it is seen from eqs. 
(3.9), (3.10). 

First there is fl quadrupole - nuclear radius interaction. It looks like 



V{x) = —rf,^Qf;WN,6^'\r). (3.11) 
Quadrupole moment of Q is related to the spin operator S of 



Q 3g(f^) 2 Q^^'^ 5 



(3.12) 

Thus contact term eq.(3.11) does not depend on the value of spin I of 
nucleus. As a result the corresponding shift of p-level can depend only on 
the value of J = L-l-S, i.e. this contact term contributes into fine splitting of 
Coulomb p-levels, not into hyper-fine ones. According to our assumption fine 
structure is mainly determined by interaction of Q magnetic moment with 
orbit. Thus interaction eq.(3.11) slightly change this fine shift of levels with 
different angular momentum j but does not split them. 

There is a similar term in interaction of nuclear quadrupole moment Qj-f ^ 
with e.m. radius of Q 
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^W = ^^(n)Q!fV.V,<5(=')(r) (3.13) 



35 'm'^ij 

with 



Qlf = 2/(2^-1) ^-^-^^- + ^^^^ - l^^^ + ^3.14) 



This interaction leads to hyperfine sphtting of p-levels indeed. The mathe- 
matics that helps to calculate matrix elements of tensor operators between 
wave functions from multiplets with different j and / is very similar to the 
case of ordinary atomic physics ^3J. The result of this calculation is following 

AE^'\nJ = = |Q(^V^^)|/(0)rX«(j)T(/,j), (3.15) 

where /(O) is the derivative of Coulomb / = 1 wave function at origin (see 
0) 

and T{f,j) and X{j) are spin-factors 



\fm' = ^^imZaf, (3.16) 



T(/, j) = IJ^{J^J^ + J^X^--j(j + l)6,^} = 2(J ■ I)2 + (J . l)_-/(/+i)j(j + l), 

(3.17) 

and 

;,»,(j)^3bb + l)- 77/12] 
83(3 + 1) 

Here 

J,/, = J ■ I = ^[/(/ + 1) - /(/ + 1) - jU + 1)]. 

For 1 = 1 nucleus and for p-states the only interesting states are those 
with j = 3/2 and j = 5/2. For this cases X^^^j = 3/2) = -4/15 and 
X(i)(j = 5/2) = 1/10. For j = 1/2, T(/, j) vanishes, as it should, smce one 
can't have a quadrupole at j = 1/2. 

In this paper we are interested mainly in Q quadrupole - nuclear 
quadrupole contact interaction 

V{x) = -^Qf Q!?[V,V, - ^V%,]5{r). (3.19) 
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Matrix elements of tensor operators for this interactions can be calculated in 
the similar way. The result is 

A(2)E(n,/ = = -lg(^)Q(^)|/(0)|2x(2)(j)T(/,j), (3.20) 

with spin factor 

\J) i6^-^^-^^ > 10 80j(j + l)^ ^ ' 

For 1 = 1 nucleus and for j9-states the only interest is in j = 3/2 and j = 5/2, 
as before. For this cases X(2)(j = 3/2) = 77/150 and X(2)(j = 5/2) = 1/50. 
For the derivation of eqs. (3.18) and (3.21) see next section. 

4 Calculation of spin-factors 

The calculation of matrix elements of tensor operators between states that 
correspond to irreducible representation of rotation group is well-known (see, 
e.g 3j). In this section we present some intermediate formulae that help to 
get eqs. (3.18), (3.21). 

Matrix elements of any tensor in the subspace of multiplets with given 
value of J can be represented as a matrix elements of appropriate combination 
of operators J. 

We start with the state with given L, namely with subspace of p-states. In 
the special basis where / = 1 wave function is represented as a vector-function 
"^{x) = x/(r) (r = |x|) one can check that 

V.V,5(=^)(r) = |/(0)|2{L25,, - L,L, - L,U} 

(z) 

Thus interaction of nuclear quadrupole moment Q]i^ with e.m. radius of 
Q can be rewritten as 

2 

Vix) = '-^QiPlfiOnU, + L,U}- (4.1) 

( Term ~ 6ij does not contribute into interaction since tensor Qlj is trace- 
less) . 

Second step is averaging of this operator over states with given J. At this 
step any tensor can be constructed from operators Jj so that 

Lit, + L,U - 2/?,LH,, = X(i)(j) ■ {Jj, + J,J, - 2/3^5,,}, (4.2) 
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where X^^\j) is an unknown constant. To calculate this constant we use the 
following relations 



To = MjJj + JjJi - -J%j}Jj = 4/3j(i + - l/2)(j + 3/2); (4.3) 
Ti = ULiLj + LjL, - h^Sij}Jj 2(J • L)2 - ( J • L) - + + 1) = 

= ^(j - l/2)(j + 3/2){j(j + 1) - 77/12}, (4.4) 

where in the last hne we take I — 1. Thus for constant X^^\j) one gets 
relation 

= Ti/To, 

that is equivalent to eq.(3.18) from the previous section. 

Quadrupole-quadrupole interaction for the states with given L can be 
rewritten as 

V{x) = -^\m\'Q\?Q\!\^L%, - LjL, - L,Lj}, (4.5) 
where the quadrupole tensors are 

= 2/(2^/ -1) ^^-^^- + - l^^^ + ^^-^^ 

and 

30(^) 2 
^^^^ ^ 2S{2S-1) ^^'^' + ~ 3^^^ + ^^^^^'^^ ^^-^^ 
For multiplet with given total angular moment J 

V{x) = -l|/(0)pgSf < q\^^{^L%, - L,Lk - L,L,} >j, (4.8) 
where < .. >j is average over multiplet with given J and 
< Qfk\^L%k - LjLk - LkLj} >j= X^^\j){jJ, + JjJ, - p%j}. (4.9) 

To calculate X^^\j) we need new relations 
T2 = Ji{SiSj+sA-p'5ij}Jj = l(j-l/2)(j+3/2){j(j + l)-7/4}, (4.10) 
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Ts — Ji{SiLj + SjLi — — S ■ L5jj}Jj — 
= - l/2)(j + 3/2){j(j + 1) + 49/4}. (4.11) 

D 

After some algebra one gets eq.(3.21) for X^'^\j). 

5 Numerical estimates 

Thus for hyperfine splitting we get 

^E{n,l = l,3j)=^-Q'~'^\fm'T{f,3){X^'\3)erl-lx^^^ (5.1) 
with 

r?2 - 1 

l/(0)^ = ^-^(mZ«)^ 

T(/ = l/2) = 5,r(/ = 3/2) = -4,T(/ 
for i = 3/2 and 

T(/ = 3/2) = 28/3, T(/ = 5/2) = -32/3, T(/ = 7/2) = 10/3 (5.4) 

for j = 5/2. One can check that hyperfine splitting satisfies well-known Sum 
Rules 

2T(l/2) + 4T(3/2) + 6T(5/2) = (5.5), 

and 

4T(3/2) + 6T(5/2) + 8T(7/2) = 0. (5.6) 

Unknown quantities in eq.(5.1) are electric radius er^ and quadrupole mo- 
ments Q'^^^ and Q^^K If we take er^ ~ 10~^/ermi^ and Q'^^^ ~ 10~^/ermz^ 
we get that for Z = 7 hyperfine splitting is of the order of keV. We believe 
that such splitting can be measured experimentally. 

The main background that interferes with hyperfine splitting is strong 
interaction of Q with nucleus ^. Little is known about strong interaction 
of Q. To get an estimate for the shift and width of Coulomb levels in Q 
atoms due to strong interaction we looked at p-levels of antiprotonic atoms 
jnj. According to these data for antiprotonic hydrogen atom T{2p) ~ lO^^eV 

^We are grateful to M.Pospelov who brought our attention to this subject. 



(5.2) 

= 5/2) = 1 (5.3) 
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and for antiprotonic deuterium atoms T{2p) --^ 5 • lO^^eV. For heavier atoms 
widths of the states grows rapidly. Thus T{2p) ~ 660eV for 6Li ^. 

We suppose that ^2 interacts with nucleus less than antiproton. Thus 
there is a good chance that strong interaction does not screen hyperfine 
structure of Q atoms due to contact interactions. It is a challenge for exper- 
imentalist to measure quadrupole moment of fi. 



6 Summary and Conclusions 

We have discussed the hyperfine spectrum of a bound system of and 
a nucleus of angular momentum one, in a scenario which corresponds to 
a small electric quadrupole moment for the of the order 10~^ fermi 
square. We find that there is a quadrupolar contact term which contributes 
to the orbital angular momentum one states. Our discussion suggests that 
contact terms are a imiversal property of interaction between multipoles, and 
not restricted only to the interaction of two dipoles- electric or magnetic. 
However, multipoles higher than quadrupoles require two particles of spin 
3/2, and the number of such systems is quite small. We assume that the 
nucleus which binds the Q~ has only an electric quadrupole moment and no 
magnetic dipole, an oversimplified case, which is approximated by the nucleus 
NitrogenlA. We do not imply that this is the optimal case to determine the 
quadrupole moment of the Q~. It remains to be found by experiment if the 
scenario we discuss is in fact useful. 
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Appendix 1. 



Electromagnetic interaction of two particles 

The interaction of two multipole moments at long distances is well known 
and can be found in text-books on classical electrodynamics. In this paper 
we are mainly interested in the contact interactions of the moments. These 
terms are less familiar. We derive them using a Feynman diagram technique 
in non-relativistic notations. 

Consider the scattering of two particles with charges (61,62), dipole mo- 
ments (di,d2), etc. The scattering amplitude T{q) is given by one photon 
exchange diagram and is equal to 

nq) = -^[p^'\q)p^'\-q) - f\q)f\-q)l (Al.l) 

where p{q) and are the matrix elements of the density and of the current 
of the particle. In the non-relativistic approximation T(g) = —U{q), where 
U{q) is the Fourier transform of the potential energy. This gives a system- 
atic way to calculate potential energy between two particles with different 
multipoles moments. 

To perform the calculation we need first of all the expression for the matrix 
element of the current operator j{q) =< P2\j\pi > for a particle with electric 
dipole moment d, with magnetic moment /x, with quadrupole moment Qij, 
etc. The expression for matrix element of electromagnetic current = (p, j) 
is the following (for derivation see Appendix 2): 

Piq) = jo{q) ^ 6[1 - ^r^q^] - idiqi - ^Qijqiqj + 
g 

ji{q) - 2(^1 + ^2)« - ieikiPiqk + (^1-2) 

where q = p2 — Pi is momentum transfer and v = p/m is velocity of a 
particle. We introduced form-factor F{q^) ~ 1 — (r^g^)/6 to take into account 
the distribution of charge. Here is the e.-m. radius of the distribution. 
The expression for the current density operator j for a particle with charge 
6 and magnetic moment /x moving in a magnetic field is well known and can 
be found in courses of quantum mechanics (see, e.g.. Landau and Lifshitz 
[3]). For the derivation of other terms in (Al.2) see Appendix 2. 
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Consider now how this procedure works in case of charge-charge inter- 
action. The scattering amphtude is 

Tiq) = -Uiq) = -^[1 - lq\rl + rl) + j-^rfry]. 

In coordinate representation l/q"^ corresponds to l/47rr, constant in g-space 
corresponds to 5-function in x-space, and q^ corresponds to the second deriva- 
tive of S- function: 



47rr ' 
l=^5(3)(r); 

te)=^-V.V/3)(r). 
As a result U {x) is given by 

Uix) = 1^ - ^-^{rl + r',)6^'\T) - ^r^^W^^)^. {A1.3) 

We get long-range Coulomb term and contact interactions of charge with 
charge radius and of two charge radii. 

Consider dipole-dipole interaction. In momentum space 

q^ 

The product q^qk can be decomposed into two irreducible tensors so that 

QiQk I QiQk 1 J \ , ^ X r Ai A\ 

= (-^ - ^Sik) + -Sik [AlA) 

Irreducible tensors in momentum space correspond to irreducible tensors 
in coordinate space. Thus Fourier transform of {qiqu/q^ — l/35jjfc) is traceless 
and symmetric second rank tensor in coordinate space made from the second 
derivatives of l/(47rr), i.e. 

- -^5ik) ^ ( ), 
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where n = r/r. Using irreducible tensors we automatically subtract con- 
tact terms ~ 6{x) from long-range potential. In spherical coordinates, the 
irreducible tensor of the second rank is proportional to the second spherical 
harmonic Y2^rn{0,4>)- 

The second term ~ 6ik in eq (A 1.4) is constant in q- space. In coordinate 
space gives 5-function, i.e. 

1^5{3)(r) 

In this way we get well known long-range potential for two electric dipoles 
in p-wave and point-like interaction of electric dipoles in s-wave 

(d.d.)-3(din)(d.n) ^l^^^^ 

Anr-^ 6 

Magnetic moment - magnetic moment interaction is very similar to 
electric moments interaction. The only difference is that /x contribute into j 
( not to p) and that ^t is axial vector so that instead of dot-products we have 
cross-products of vectors: 

Uiq) = — \{eikiqkP^j^^^){evmnqmPn^) (Al.6) 
Transforming this expression to coordinate space we find 

^^^^ ^ /^S'Vr-3(/iPnO(/xi^^n.) _ 2^^(i)^(2))^(3)(^) 
Aur-^ 3 

Long-range interaction is known from magneto-statics. The result for contact 
term was found originally by Fermi and independently by Casimir in the 
framework of relativistic Dirac equation. Later it was derived using non- 
relativistic Breit equation in the book [Zj. Elementary derivation of contact 
terms in magneto-statics is also known (see [S] ). 

We are now ready to derive contact terms for quadrupole interaction. 

Quadrupole-electric radius interaction in momentum space looks as 

V{q) = ^[eirfQlf + e2rlQff]q.q,. {A1.8) 
In coordinate space it gives 

V^(^) = -^[e,rlQ\f + e2rlQ\f]V.Vj5^'\v). (A1.9) 
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This interaction vanishes for s-wave but is nonzero for p-wave. 
Quadrupole-quadrupole interaction in momentum space has the fol- 
lowing form 

= (ALIO) 

where we have to express the product of the four momentum factors in terms 
of irreducible tensors. Different tensors represent interaction in d-, p- and 
s-wave. 

We construct irreducible tensor Tj^jli as a sum of all possible symmetric 
terms: 

T-^li = qiqjqkqi + M^ikjquqi + kiqiqj + kkqjqi + Sjiqiqu + Suq^qk + Sjkqiqi)+ 

+Bq'^(5ij5ki + 5ik5ji + 5ii5jk)- 

Irreducible tensor T^jli has to be traceless. The constraint T^^,^i = gives 
A^-{l/7) and B = (1/35). 

Thus eq. (A 1.10) can be rewritten in following form 

In coordinate space irreducible tensor Tj^jli is proportional to fourth harmonic 
^4,m(^5 0) • Thus this term in eq. (A 1.11) represents the well-known long- 
range interaction of two quadrupoles in d-wave (see eq. (2.1)). 
When we go to coordinate space 

iW) - 



47rr 



{qiqj)^-ViVjS^'\r). 

In this way differentiating 1/r we get the standard formula for quadrupole- 
quadrupole interaction at long distance. Second term in cq. (A 1.11) is 
quadratic polynomial in momentum space. In coordinate space it corresponds 
to second derivative of delta function 

Vcontix) = -^gg^Q;?{V,V, - ^V%k}S^'\r), {A1.12) 

which is equivalent to eq.(2.2) used in the text. This result is new. 
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Appendix 2. 



Electromagnetic vertex 

One can find the standard derivation for tlie current density operator j 
for a particle with charge e and magnetic moment /x moving in a magnetic 
field in quantum mechanics text-books (see, e.g., Landau and Lifshitz [3]). 
The matrix element of this operator between two states with momenta pi 
and p2 is 

" 6 

< P2\ji\Pi >= -{vi + V2)i + ieikiQkl^h (^2.1) 

where q = p2 — Pi is momentum transfer and v = p/m is velocity of a 
particle. For the charge density p = jo the matrix element is 

< P2\p\pi >= eF{q^) ^ e(l - -r^g^). (A2.2) 

6 

We need similar expression for the electromagnetic current of a particle 
with electric dipole moment d, with quadruple moment Qij ( and with mag- 
netic moment as well). One can derive this expression classically without 
any reference to quantum mechanics. 

That is two-steps procedure (see also ref. [IHl)- First, notice that the 
interaction of the particle with electromagnetic field = (0, A) is described 
by currents 

H,nt = P(l) - iA. (A2.3) 

This is the definition of the current. 

Second, notice that in the static field the interaction energy of a particle 
is given by the multipole expansion. By definition of the multipole moments 
d, n, Qij, etc, we have 

1 OE ■ 

Hint ^ e0(O) - diEi - -Qij^ - PiHi = 

6 oxi 

= [e + diVi + ^QijViVj]</)(x)|a;=o - eikiPiV kAi{x)\^=o, {A2A) 

where Ei{x) = — Vi0(a;), Hi = CikiV kAi{x) are electric and magnetic fields 
respectively. 

Comparing eqs.(A2.3) and (A2.4) we find the contribution of electric 
dipole, magnetic dipole and quadrupole moments into e.-m. current. 



17 



In coordinate space it is equal to the sum of differential operators 
jo = e + diVj + ^QijVjVj + 

ji = eikiHiVk- (^2.5) 

In momentum space 

jX?) = J j{x)exp{—iqx)d^x 

it can be written as 

jo = e - idiQi - ^Qijqiqj + 

ji = -ieikiHiqk + ■■■■ {A2.6) 

Fourier transform of the classical current corresponds to the matrix ele- 
ment of the quantum operator j(g) =< P2\j\pi >■ 

Thus we get the complete representation of the matrix element of the 
current operator in momentum space 

jo{q) ~ e[l - ^r^g^] - idiqi - ^Qijqiqj + 

g 

JiiQ) - 2 + ^2); - ieikil^iqk + ■■■■ {A2.7) 

In the language of Feynman diagrams eq.(A2.7) corresponds to the vertex 
operator that describes interaction of the particle with photon. 
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